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Recent theory predicted that the Quantum Spin Hall Effect, a fundamen- 
tally novel quantum state of matter that exists at zero external magnetic field, 
may be realized in HgTe/(Hg,Cd)Te quantum wells. We have fabricated such 
sample structures with low density and high mobility in which we can tune, 
through an external gate voltage, the carrier conduction from n-type to the 
p-type, passing through an insulating regime. For thin quantum wells with 
well width d < 6.3 nm, the insulating regime shows the conventional behavior 
of vanishingly small conductance at low temperature. However, for thicker 
quantum wells (d > 6.3 nm), the nominally insulating regime shows a plateau 
of residual conductance close to 2e^/h. The residual conductance is indepen- 
dent of the sample width, indicating that it is caused by edge states. Further- 
more, the residual conductance is destroyed by a small external magnetic field. 
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The quantum phase transition at the critical thickness, d — 6.3 nm, is also in- 
dependently determined from the magnetic field induced insulator to metal 
transition. These observations provide experimental evidence of the quantum 
spin Hall effect. 
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The theoretical prediction of the intrinsic spin Hall effect in metals and insulators (|i]|2] 

has generated great interests in the field of spintronics, since this effect allows for direct 
electric manipulation of the spin degrees of freedom without a magnetic field, and the resulting 
spin current can flow without dissipation. These properties could lead to promising spintronic 
devices with low power dissipation. 

However, beyond the potential technological applications, the intrinsic spin Hall effect has 
guided us in the search for new and topologically non-trivial states of matter. The quantum Hall 
state gives the first, and so far the only example of a topologically non-trivial state of matter, 
where the quantization of the Hall conductance is protected by a topological invariant. The 
quantum spin Hall (QSH) insulators (|4]|5]|^ have a similar, but distinct non-trivial topological 
property. The QSH insulators are invariant under time reversal, have a charge excitation gap in 
the bulk, but have topologically protected gapless edge states that lie inside the bulk insulating 
gap. This type of insulator is typically realized in spin-orbit coupled systems; the corresponding 
edge states have a distinct helical property: two states with opposite spin-polarization counter- 
propagate at a given edge (|?][^|^. The edge states come in Kramers' doublets, and time reversal 
symmetry ensures the crossing of their energy levels at special points in the Brillouin zone. 
Because of this energy level crossing, the spectrum of a QSH insulator cannot be adiabatically 
deformed into that of a topologically trivial insulator without helical edge states; therefore, in 
this precise sense, the QSH insulators represent a topologically distinct new state of matter. 

It has been proposed theoretically that HgTe/(Hg,Cd)Te quantum wells provide a natural 
realization of the quantum spin Hall effect (6). In zincblende-type semiconductor quantum 
wells, there are four relevant bands close to the Fermi level. The El band basically consists of 
the two spin states of the s orbital, while the HHl band basically consists of the \px + ipy, T) 
and I — (p,j. — ipy), I) orbitals. The effective Hamiltonian near the T point is given by 

Heffih, ky) = (^^^^ ) , if = eik) + d,{k)a, (1) 
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where are the Pauli matrices, and 

di + id2 = A{kx + iky) = Ak+ 
d3 = M- B{kl + kl),ek = C- D{kl + kl). (2) 

Here, k^ and ky are momenta in the plane of the 2DEG, while A, B,C, D are material specific 
constants. Spin-orbit coupling is naturally built-in in this Hamiltonian through the spin-orbit 
coupled p orbitals \px + ipy,1) ^^d \ — {px — ipy), i)- Two dimensional materials can be 
grouped into three types according to the sign of the Dirac mass parameter M. In conventional 
semiconductors such as GaAs and CdTe, the s-like El band lies above the p-like HHl band, 
and the mass parameter M is positive. Semi-metals such as graphene (|9 170) are described by a 



massless Dirac model with M = 0, although the bands have a different physical interpretation. 
In so-called "inverted" semiconductors such as HgTe, the s-like orbital lies below the p-like 
orbitals; therefore, the Dirac mass parameter M in the HgTe/(Hg,Cd)Te quantum wells can be 
continuously tuned from a positive value M > for thin quantum wells with thickness d < d^ 
to a negative value Af < for thick quantum wells with d > d^.. A topological quantum phase 
transition occurs atd = dc, where the system is effectively described by a massless Dirac theory 
just like for graphene. The nature of this quantum phase transition has also been investigated 
in more realistic models beyond the simple four band model presented here, reaching the same 
conclusion ( [77] ). 

The QSH phase occurs in the inverted regime where M < 0, i.e., when d > dc- The sample 
edge can be viewed as a domain wall of the mass parameter M, separating the topologically 
non-trivial phase with M < from the topologically trivial phase with M > 0, which is adia- 
batically connected to the vacuum ([72]). Massless helical states are confined on the sample edge. 
The sample has a finite conductance even when the Fermi level lies inside the bulk insulating 
gap. Therefore, as suggested in Ref. [6] the QSH state can be experimentally detected by mea- 



suring a residual conductance plateau as one varies the gate voltage in the nominally insulating 
regime. Furthermore, because the current is carried by the edge states, the conductance should 
be independent of sample width. Protected by the time reversal symmetry, non-magnetic impu- 
rities or any other time-reversal invariant local perturbations cannot cause elastic backscattering 
of the helical edge states, which warrants the topological robustness of the edge state conduc- 
tance. However, the presence of magnetic field breaks time reversal symmetry, therefore can 
open up a gap in the energy spectrum of the edge states, and remove the residual conductance 
due to the edge states. 

We set out to test these theoretical predictions by measuring the transport properties of 
HgTe/(Hg,Cd)Te quantum wells as a function of the sample thickness, the gate voltage and the 
external magnetic field. We use modulation-doped type-III ( [TJ] ) HgTe/Hgo.sCdo.yTe quantum 



well (QW) structures fabricated by molecular beam epitaxy with widths p5) varying 
from 5.5 nm (d < dc) to 12 nm {d > d^). Dedicated low-thermal budget optical and e-beam 
lithography is used to structure devices in Hall bar geometry with dimensions, (L x W), of 
(600 X 200) /im^, (20.0 x 13.3) /im^, (1.0 x 1.0) /xm^, and (1.0 x 0.5) /im^ (see inset in Fig. 



lA). All devices have a 110 nm thick Si3N4/Si02 multilayer gate insulator (16) and a 5/30 nm 
Ti/Au gate electrode stack. Transport measurements are done in a '^He/^He-dilution refrigerator 
(base temperature T < 30 mK, uniaxial fields up to 18 T) and in a ^He cryostat fitted with 
a vector magnet (T = 1.4 K, and fields up to 300 mT in variable direction), using lock-in 
techniques. At zero gate voltage, the samples are ra-type, exhibiting carrier densities between 
1.3 X 10^^ cm"^ and 3.5 x 10^^ cm"^ and mobilities up to 1.5 x 10^ cm^/(Vs). The carrier 
density can be reduced continuously by applying a negative gate voltage to the Au electrode 
with respect to carriers in the QW. The Si-O-N gate insulator stack allows for quite large gate 
voltages, enabling us to gate the samples through the gap from n-type to p-type conductance. 
The change in carrier type can be monitored from Hall experiments, as shown in Fig. lA 
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for a large [{L x W) = (600 x 200) /im^] Hall bar with a well width of 6.5 nm, at 30 mK. The 
change in carrier type is directly reflected in a sign change of the slope of the Hall resistance 
Rxy, and we can directly infer that the carrier density varies from n = 1.2 x 10^^ cm^^ at gate 
voltage =-1.0Vtop = 1.0 x 10^° cm"^ at = -2.0 V. At modest magnetic fields, for both 
n-type and p-type channels, R^y exhibits quantum Hall plateaus, indicative of the good quality 
of the material, until at fields above ca. 3 T the last Landau level is depleted. 

Remarkable transport behavior is observed for -1.9 V < Vg < -1.4 V (see in particular the 
green and the red traces in Fig. 1), where the sample is insulating at zero magnetic field (i.e., 
the Fermi level is in the gap). For these gate voltages, we observe that the sample undergoes 
a phase transition from an insulating state to a QH state with a quantized Hall conductance of 
Gxy = ±e^//i, either n- or p-type depending on Vg, at a small (order 1-2 T) applied magnetic 
field. The sample remains in the QH state for a few more T, and then becomes once again 
insulating. We have observed this phenomenon in a number of samples in the inverted regime, 
with 6.5 nm < d < 12 nm. 

The phase transition from an insulating state to a QH state is a non-trivial consequence of the 
inverted band structure of the QSH insulator, and can be explained by the level crossing of the 
El and HHl Landau levels, which can be directly obtained from the minimal coupling of the 
simple Dirac model Eq. [T]to a perpendicular magnetic field B±. If we only consider the orbital 
effects of the magnetic field, two series of Landau levels are obtained from the upper and lower 
2x2 blocks of Hamiltonian, Eq.[T} in which the two levels closest to the Fermi energy are given 
hyE+ = C + M-{D + B)l-^ and E_ = C - M - {D - B)l-\ with = ^h/{eB^). Thus 
the condition for level crossing is given by E+-E_ = 2M ~ 2Bl-^ = or 5^ = {hM) /{eB) . 
Generally, the B parameter is always negative, therefore, we can see that the level crossing 
occurs only in the inverted region with M < 0. 

The Landau levels for the normal {d < dc) and the inverted {d > dc) regime are shown 
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in Fig. 2A and B, respectively. Edge states in the presence of an external magnetic field can 
be easily obtained by solving our simple Dirac model in Eq. [T] with open boundary condition 
along one direction and periodic boundary condition along the other direction. Fig. 2C shows 
the bulk and edge states for a conventional insulator. With increasing thickness d, the two states 
closest to the Fermi energy approach and then cross each other. This "band inversion" leads 
to the bulk and edge states [Fig. 2D]. The Fermi energy crosses a pair of counter-propagating 
edge states on each edge, resulting in no net Hall effect. These counter-propagating edge states 
are similar to the helical edge states of the QSH insulator. However, due to the presence of the 
magnetic field and the breaking of the time reversal symmetry, they are not robustly protected. 
Starting from this case, increasing magnetic field shifts the red Landau level towards higher 
and the blue one towards lower energies. When one of the bulk Landau level crosses the Fermi 
level, there is a single edge state on each edge left, and one obtains a net QH effect with either 
Gxy = e^/h [Fig. 2E], or with G^y = —e^/h [Fig. 2F]. When the magnetic field is increased 
further, the second bulk Landau level crosses the Fermi level, and one reaches the conventional 
insulator case [Fig. 2C] but with the colors of the Landau level interchanged. In models with 



bulk inversion asymmetry (BIA) {11 1, the level crossing between El and HHl Landau levels at 
Sj_ can be avoided, and the phase regions (i) and (ii) in Fig. 2B become connected. Generally, 
the non-monotonic dependence of the Landau level energies leads to the transition from the 
insulating state to the QH state at a constant Fermi energy, when the magnetic field is varied. 

While the four band Dirac model (Eq. [T]) gives a simple qualitative understanding of this 
novel phase transition, we have also performed more realistic and self-consistent eight band 



k ■ p model calculations (13) for a 6.5 nm quantum well, with the fan chart of the Landau levels 
displayed in Fig. IB. The two anomalous Landau levels cross at a critical magnetic field B^, 
which evidently depends on well width. This implies that when a sample has its Fermi energy 
in the gap at zero magnetic field, this energy will always be crossed by the two anomalous 
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Landau levels, resulting in a QH plateau in between the two crossing fields. Fig. 3 summarizes 
the dependence of 5^ on well width d. The open red squares are experimental data points that 
result from fitting the 8-band k • p model to experimental data as in Fig. 1, while the filled red 
triangles solely result from the k ■ p calculation. For reference, the calculated gap energies are 
also plotted in this graph as open blue circles. The band inversion is reflected in the sign change 
of the gap. Note that for relatively wide wells (d > 8.5 nm) the (inverted) gap starts to decrease 
in magnitude. This is because for these well widths, the band gap no longer occurs between 
the El and HHl levels, but rather between HHl and HH2 - the second confined hole-like 



level, as schematically shown in the inset of Fig. 3 [see also (17)]. Also in this regime, a band 
crossing of conductance- (HHl) and valence- (HH2) band derived Landau levels occurs with 



increasing magnetic field ( |i5][i^[iq ). Fig. 3 clearly illustrates the quantum phase transition that 
occurs as a function of d in the HgTe QWs: only for d > dc, exists, and at the same time 
the energy gap is negative (i.e., the band structure is inverted). The experimental data allow for 
a quite accurate determination of the critical thickness, yielding d^ = 6.3 ±0.1 nm. 

The actual existence of edge channels in insulating inverted QWs is only revealed when 
studying smaller Hall bars (note that the typical mobility of 10^ cm^/(Vs) in n-type material 



implies an elastic mean free path of /mfp ~ 1 /um (19 20) - and one may anticipate lower 
mobihties in the nominally insulating regime). The pertinent data is shown in Fig. 4, which 
plots the zero 5-field four terminal resistance -Ri4,23 = V23//14 as a function of normalized gate 
voltage (Ythr is defined as the voltage for which the resistance is largest) for several devices, that 
are representative for the large number of structures we have investigated. -Ri4,23 is measured 
while the Fermi level in the device is scanned through the gap. In the low resistance regions at 
positive Vg — Vthr the sample is n-type, at negative Vg — Vthr the sample is p-type. 

The black curve labeled I in Fig. 4 is obtained from a medium sized [(20.0 x 13.3) /im^] 
device with 5.5 nm QW, and shows the behavior we observe for all devices with a normal band 
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structure: when the Fermi level is in the gap, -Ri4,23 increases strongly and is at least several tens 
of Mfi (this actually is the noise floor of the lock- in equipment used in the experiment). This 
clearly is the expected behavior for a conventional insulator. However, for all devices containing 
an inverted QW, the resistance in the insulating regime remains finite. -Ri4,23 plateaus out at 
well below 100 kfi (i.e., 6*14,23 = 0.3e^//i) for the blue curve labeled II, which is again for a 
(20.0 X 13.3) fxm^ device fabricated by optical lithography, but now contains a 7.3 nm wide 
QW. For much shorter samples (L = 1.0 /im, green and red curves III and IV) fabricated from 
the same wafer, ^14,23 actually reaches the predicted value close to 2e^ /h. This observation 
provides firm evidence for the existence of the quantum spin Hall insulator state for inverted 
HgTe QW structures. 

In Fig. 4, we have included data on two c? = 7.3 nm, L = 1.0 /im devices. The green trace 
(III) is from a device with W = 1.0 /im, the red trace (IV) corresponds to ly = 0.5 /xm. Clearly, 
the residual resistance of the devices does not depend on the width of the structure, which is 



evidence that the transport occurs through edge channels (21 ). One notices that the traces for 
the d = 7.3 nm, L = 1.0 yum devices do not reach all the way into the p-region. This is because 
the electron-beam lithography needed to fabricate the devices increases the intrinsic (Vg = 
V) carrier concentration. In addition, one notices fluctuations on the conductance plateaus in 
traces II, III, and IV. These fluctuations are actually reproducible and do not stem from, e.g., 
electrical noise. While all -Ri4,23 traces discussed so far were taken at base temperature (30 
mK) of our dilution refrigerator, the conductance plateaus are not at all limited to this very low 
temperature regime. In the inset of Fig. 4 we reproduce the green 30 mK trace III on a linear 
scale, and compare it with a trace (in black) taken at 1.8 K from another (L x W) = (1.0 x 
1.0) /im^ sample, that was fabricated from the same wafer. In the fabrication of this sample, 
we used a lower illumination dose in the e-beam lithography, resulting in a better (but still not 
quite complete) coverage of the n-i-p transition. Clearly, in this further sample, and at 1.8 
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K, the 2e^//i conductance plateau is again present, including (thermally smeared) conductance 
fluctuations. 

In the pure two terminal geometry, with only source and drain contacts (contacts 1 and 
4, inset Fig lA), the two counter propagating helical edge states at one given edge connect 
the chemical potential from the source and drain, respectively, and they are not in equilibrium 
with each other because the elastic backscattering vanishes between these two channels. In 
the absence of voltage probes 2,3,5, and 6, as indicated in the inset of Fig. 1, the two terminal 
conductance should give 2e^ /h. Now we consider the presence of the voltage probes. Since 
our voltage probes are not spin sensitive, the voltage measurement necessarily leads to the 
equilibration of the two helical channels with the opposite spin orientation. A simple Landauer- 
Biittiker type of calculation shows that the four terminal resistance should in fact be given by 
-Ri4,23 = h/2e^. In the presence of the voltage probes, the voltage drops V12, V23 and V34 add 
in series to give a higher resistance of Ru = V14/J14 = 3h/2e^. These results are valid as long 
as the distance between the voltage probes is less than the inelastic mean free path I in. While 
elastic scatterers can not cause backscattering of the helical edge states, inelastic scatterers can. 
We estimate the inelastic mean free path to be /j„ > l/im at our measurement temperature. 
Therefore, for the large sample (II), where the distance between the voltage probes exceeds the 
inelastic mean free path /j„, we expect the residual resistance to be higher, consistent with the 
experimental measurement shown in the trace (11) in Fig. 4. 

Another intriguing observation is that the QSH effect is destroyed by applying only a small 
magnetic field perpendicular to the 2DEG plane. Fig. 5 shows that the magnetoresistance is 
actually strongly anisotropic. (This data has been obtained in the vector magnet system at 1.4 
K.) A very sharp, cusp-like conductance peak is observed for perpendicular field, with the full 
width half-maximum (FWHM) B^^^'^ ~ 28 mT (22). The peak broadens strongly when the 
magnetic field is tilted into the QW plane. For fully in-plane fields, the QSH conductance can 
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be observed over a much wider magnetic field range (iJ^WRM ^ 0.7 T). 

The robustness of the helical edge states is ensured by the time -reversal symmetry. A mag- 
netic field breaks time reversal symmetry, and thus turns on a gap between the two otherwise 
degenerate helical edge states. The perpendicular and in-plane magnetic field lead to different 
gaps, depending on the matrix elements of the magnetization operator: 

^gap± = (Tl (z ■ f X } + g±iJ,BS_i_^ li) \B\ and Eg^p\\ = {]\ g\\iJ,BS\\ |i) \B\, in which r,j are the 
position and electric current operator, respectively, and z is the unit- vector perpendicular to the 
quantum well plane. S±qi-) stands for the dimensionless part of the Zeeman-coupling matrix 
element in perpendicular (parallel) magnetic field. We can estimate the magnitude of these two 
gaps by noting that (t| z • f x j ||) ~ ev^ and (t| S\\^± ||) ~ 1, in which v and ^ are the Fermi 
velocity and width of the edge channels, respectively, v and ^ can be obtained from the Dirac 
parameters as v = A/h and ^ ^ hv/ \ M\. The parameters for the d = 7.3 nm quantum well 
give the dimensionless ratio ev^/fiB ~ 280, which thus leads to -E'gap±/-E'gap|| ~ 10^. From 
this estimate, we can see that the strong anisotropy observed in the experiments originates from 
the high Fermi velocity of the edge states and the small bulk gap M, which together make the 
orbital magnetization dominant. 

So far our experiments have only measured the charge transport properties. While the QSH 
effect manifests itself already in the change in transport properties, we would still like to ex- 



perimentally confirm the spin accumulation resulting from the spin Hall effect (23 24) in the 
topologically non-trivial insulating regime, and compare both electric and magnetic results with 
the experiments of the spin Hall effect in the metallic regime. It would also be interesting to 
explore the regime close to the quantum phase transition point oi d = d^, and compare the 
transport properties with that of the recently discovered graphene. In many ways, the HgTe 
quantum well system can be viewed as a tunable graphene system, where the Dirac mass term 
can be tuned continuously to zero from either the positive (topologically trivial) or the negative 
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(topologically non-trivial) side. 
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1 Figure Captions 

Fig. 1 (A) Hall resistance, R^y, of a (L x ly) = (600 x 200) /im^ QW structure with 6.5 nm 
well width for different carrier concentrations obtained for gate voltages of -1.0 V (black), 
-1.1 V (purple), -1.2 V (navy), -1.3 V (blue), -1.35 V (cyan), -1.4 V (green), -1.7 V (red), 
-1.8 V (orange), -1.9 V (brown), and -2.0 V (black, lower curve). For decreasing Vg, the n- 
type carrier concentration decreases and a transition to a p-type conductor is observed, passing 
through an insulating regime between -1.4 and -1.9 V at i? = T. The inset shows a schematic 
sample layout with ohmic contacts labeled 1 to 6. The grey shaded region indicates the top gate 
electrode and defines the part of the sample where the carrier concentration and type can be 
changed. Red and blue arrows indicate the counter-propagating spin-polarized edge channels 
of the QSH effect. (B) The Landau level fan chart of a 6.5 nm quantum well obtained from 
an eight band k • p calculation. Black dashed lines indicate the energetic position of the Fermi 
energy, Ep, for Vg = —1.0 and —2.0 V. Red and green dashed lines correspond to position 
of the Fermi energies of the red and green Hall resistance traces of Fig. A. The Landau level 
crossing points are marked by arrows of the same color. 

Fig. 2 Bulk and edge state spectrum of the four band Dirac model described by Eq. [T] in 
the presence of an external orbital magnetic field. (A) The bulk Landau levels in the normal 
regime. (B) The bulk Landau levels in the inverted regime. A pair of low lying Landau levels 
cross at a finite magnetic field B^. The crossing divides the phase diagram of gate voltage and 
magnetic field into four regimes, labelled (i-iv) in the figure. (C) The low lying bulk and edge 
state energies as a function of the centers of the Landau orbitals in the normal regime. (D) The 
low lying bulk and edge state energies as a function of the centers of the Landau orbitals for the 
inverted regime, where the Fermi energy lies in between the two bulk inverted Landau levels. 
The Fermi energy crosses the Landau levels, giving rise to the one pair of counter-propagating 
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edge states. When the magnetic field is increased, the two lowest lying bulk Landau levels 
approach each other, and they cross the Fermi energy in different order, depending on the value 
of the gate voltage. The crossing of the Fermi energy by one of the Landau levels gives rise to 
either the n— or the p-type QHE for the cases shown in figures (E) and (F). 

Fig. 3 Crossing field, B^, (red triangles) and energy gap. Eg, (blue open dots) as a function 
of QW width d resulting from an eight band k ■ p calculation. For well widths larger than 
6.3 nm the QW is inverted and a mid-gap crossing of Landau levels deriving from the HHl 
conductance and El valence band occurs at finite magnetic fields. The experimentally observed 
crossing points are indicated by open red squares. The inset shows the energetic ordering of the 



QW subband structure as a function of QW width d. [taken from Ref. (17)] 



Fig. 4 The longitudinal four-terminal resistance, -Ri4,23, of various normal {d = 5.5 nm) 
(I) and inverted (d = 7.3 nm) (II, III, and IV) QW structures as a function of the gate voltage 
measured for i? = T at T = 30 mK. The device sizes are (20.0 x 13.3) fxm^ for device I and 
II, (1.0 X 1.0) yum^ for device III, and (1.0 x 0.5) jjm^ for device IV. The inset shows -Ri4,23(^g) 
of two samples from the same wafer, having the same device size (III) at 30 mK (green) and 1.8 
K (black) on a linear scale. 

Fig. 5 Four-terminal magnetoconductance, 6*14,23, in the QSH regime as a function of tilt 
angle between the plane of the 2DEG and applied magnetic field for ad = 7.3 nm QW structure 
with dimensions (L x W) = (20 x 13.3) fim^ measured in a vector field cryostat at 1.4 K. 
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